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Spectrum analyzer having a resolution filter that can be 
adjusted via a phase- variation parameter 

FIELD OF THE INVENTION 

The invention relates to a resolution filter for a 
spectrum analyser. 

BACKGROUND OF THE INVENTION 

In spectrum analysis, a specified frequency range is 
"swept" with a resolution filter with a specified 
bandwidth. The resolution filter is therefore also 
sometimes referred to as a sweep filter. A resolution 
filter of this kind for a spectrum analyser with an 
analogous structure is known, for example, from US 
5,736,845. With resolution filters of the known, 
analogous structure, only a limited sweep rate can be 
achieved; the so-called K-factor, which indicates the 
rate of sweep, is limited in resolution filters of the 
known structure. 

It has been assumed hitherto in the context of spectrum 
analysis that it is possible to sweep within a T res in the 
order of magnitude around 1/B res = T res , so that the 
resolution filter can still build up. However, this 
statement is only correct if a fixed filter is assumed 
for all the sweep rates. 

A digital resolution filter for a spectrum analyser is 
known from DE 101 05 258 Al . The resolution filter 
described in this context is characterised by a Gaussian 
impulse response. It is a so-called linear-phase 
resolution filter. Linear-phase resolution filters have a 
relatively long group delay response. As a result, during 
sweeping, these filters have a considerable frequency 
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overshoot, and the middle of the spectrum is no longer 
disposed at the frequency origin. A degree of design 
freedom, which would allow a compensation of these 
undesirable effects is not provided with the impulse 
response of the resolution filter defined in DE 101 05 
258 Al. 

SUMMARY OF THE INVENTION 

A need therefore exists for providing a spectrum analyser 
and an associated resolution filter, wherein the impulse 
response of the resolution filter has a free design 
parameter, which allows a compensation of the frequency 
overshoot, the displacement of the frequency origin and 
other undesirable effects. 

According to an embodiment of the invention, the free 
variation parameter k 0 is introduced into the phase factor 
of the impulse response. This free variation parameter 
represents a degree of freedom of the phase in the design 
of the filter. Accordingly, for example, not only linear- 
phase, but also minimal-phase filters can be realised in 
an efficient manner. 

The free variation parameter ko can preferably be set in 
such a manner that the frequency overshoot determined by 
the group delay of the resolution filter is compensated. 

As an alternative or at the same time, the variation 
parameter ko can also be set in such a manner, that the 
middle of the frequency response of the resolution filter 
is disposed at the frequency origin, that is to say, at 
the frequency f =0 . 
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Still other aspects, features, and advantages of the 
present invention are readily apparent from the following 
detailed description, simply by illustrating a number of 
particular embodiments and implementations, including the 
best mode contemplated for carrying out the present 
invention. The present invention is also capable of 
other and different embodiments, and its several details 
can be modified in various obvious respects, all without 
departing from the spirit and scope of the present 
invention. Accordingly, the drawing and description are 
to be regarded as illustrative in nature, and not as 
restrictive . 

BRIEF DESCRIPTION OF THE DRAWINGS 

The invention will be explained in greater detail below 
with reference to the drawings. The drawings are as 
follows : 

Figure 1 shows a block circuit diagram of a spectrum 



analyser, in which the resolution filter 
according to the invention can be used; 



Figure 2 



shows a block circuit diagram of the spectrum 
analysis in the equivalent baseband; 



Figure 3 



shows the impulse response of a linear-phase 
filter and of a minimal-phase Gaussian filters- 



Figure 4 



shows a sweep with a minimal-phase filter and 



Figure 5 



shows a sweep with a linear-phase filter. 



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
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Figure 1 provides an overview of a spectrum analyser 20, 
in which the resolution filter 29 according to the 
invention is used. Figure 1 shows only the region of the 
signal below the intermediate frequency level, which is 
relevant here. 

The intermediate frequency signal marked ZF in the 
drawings is filtered in a band-pass filter 21. The band- 
pass filter 21 is connected to an analog/digital 
converter 22. This is followed by the I/Q mixing 23 in an 
I/Q demodulator 24, which, in the conventional manner, 
comprises a local oscillator 25 with two outputs phase- 
displaced through 90°, which are supplied, together with 
the filtered and A/D-converted intermediate frequency 
signals, to a mixer 27 of the I-branch and a mixer 26 of 
the Q-branch respectively. 

This stage is followed by a digital filtering 28 with the 
resolution filter 29 according to the invention. An 
envelope-curve ic^Lif ication 21 then takes pla.np in an 
envelope-curve rectifier 32. Logging 33 is implemented in 
a logging unit 34. The logging unit 34 is connected to a 
video filter 36, in which a video filtering 35 is 
implemented . 

Various detectors 38 to 41, for example, a peak detector 
38, an auto peak detector 39, a sample detector 40 and an 
RMS (Route Mean Square) detector may be provided in the 
detection stage 37. Depending on the requirements, either 
all four detectors can be incorporated in a high- 
performance spectrum analyser 20, or only specific 
detectors, for example, only a single detector for 
specialised measuring tasks, needs to be installed. 
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Evaluation and control are implemented via a 
microprocessor 42 . 

Figure 2 shows a simplified block circuit diagram of 
blocks 24, 29 and 32 of the spectrum analyser 20 from 
Figure 1. The complex input signal v(t) to be analysed is 
supplied to a complex-conjugate former 2, which forms the 
complex-conjugate signal v* (t) of the input signal v(t). 
The complex-conjugate input signal v*(t) is mixed down 
into the baseband signal x(t) in a mixer 3 by 
multiplication with the sweep signal e j(p<f> . The frequency 
f(t) of the sweep signal is shown at the top of Figure 2 
as a function of time t, wherein it can be seen, that the 
sweep frequency f (t) changes in a linear manner with 
time. The phase angle <p(t) as a function of time t is 
obtained by integration. The baseband signal x(t) is 
supplied to the resolution filter 4 according to the 
invention. In the resolution filter 4, the baseband 
signal x(t) is convoluted with the impulse response 
hused(t) of the resolution filter 4. This leads to the 
output signal y(t). The modulus iy(t) | of the signal y(t) 
is formed in a modulus former 5. 

The lower region of Figure 2 shows by way of example an 
input signal v(t), of which the spectrum consists of two 
discrete spectral lines. An example of the transmission 
function H(t) of the resolution filter 4 is also shown. 
The spectrum shown to the right of this is displayed at 
the output of the spectrum analyser 1, wherein the 
spectral lines are widened by the resolution bandwidth 
B res of the resolution filter 4. The resolution bandwidth 
B res corresponds to the bandwidth with an attenuation by - 
3dB relative to the maximum. 
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By way of further explanation of the invention, the 
deliberations from document DE 101 05 258 Al, which lead 
to a resolution filter with a given impulse response, are 
briefly presented again below. 

The spectrum of the signal v(t) is initially windowed 
with the impulse response of the resolution filter, and 
following this, the Fourier transform is implemented 
according to the following equation: 



S(f) = J v(T)h res (T)-e- je)T dr = H res (f) * V(f) 

-00 

(1) . 

The question of the correlation of the spectrum with 
white noise is relevant in this context. This correlation 
describes the distance, at which the spectrum is un- 
corrected. The ACF (auto-correlation function) of the 
input signal with white noise is described by the 
following equation : 

E{v(r)v*(r + % -NJ2 S(dt) 

real/imag 

(2) . 

The ACF of the Fourier spectrum is obtained using 
equation (1) . 

{ 00 00 
J v (r,)^,) dx x ■ J v (r 2 )h res (r 2 )-e- J ^ d ^ dr 2 
-00 -00 

= J ] ^(rO-v^^fr^.^Cr^e-^-^V^^ dr x dr 2 



7 



Substitution of equation (2) with Ti=t 2 :=t leads to the 
following equations: 

E{S\/) • S(f + df)} = J 7V 0 ti res {r) • /z rej (r) e"^ r 

-00 

-00 

= N,-F{\h res {rf) 

For a Gaussian filter, the following applies: 



h (t) = \^—B >e 

gauss V ) -y21n(2) re5 



21n(2) 



r^ 2 



\ ! res J 



h ( n = 

~~ gauss w s 



-21n(2> 



P. 



B 



res J 



(3) 



The equations below follow from equation (3) : 



F- X {\H gauss {fi\ 



= F~\e 

= i 



-21n(2>2- 



-B 

21n(2) res 



21n(2) 



= P^-B„Jyl2-e 

\21n(2) 



21n(2) 



with5 w = 5 w /V2 
with C - 7/^ 



(4) 



8 



Moreover, with equation (4), the equations below follow: 




With a Gaussian resolution filter, the equation below is 
obtained using equation (5) : 



-Inf2vW 2 

E{s\f)-S(f + df)\ = N 0 -B mlse -e 

(6) . 

Figure 2 shows a block circuit diagram of a spectrum 
analysis in the equivalent baseband. It should be borne 
in mind that the HF-signal v(t) to be investigated is 
observed in the equivalent baseband in order to allow a 
simpler model (i.e. without spectral components at f<0). 
After the formation of v*(t), the signal is multiplied by 
the rotary phasor e j(p(t) leading to the following equation: 



x(t) =v(t)-e Mt) 



(7) 
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The frequency of the rotary phasor increases in a linear 
manner with time according to the following equation: 



The K-factor indicates the sweep rate. Since the 
resolution filter requires a transient time of 
approximately T res , the frequency should change by a 
maximum of B res within T reS f which, according to equation 
(8), corresponds to a maximum k-factor of k=l. The 
following phase is obtained by integration: 



The signal x ( t ) is then f i 1 1-ptpH through the resolution 
filter with the impulse response h use d(t) producing the 
output signal y(t). The envelope curve ly(t)| is 
determined from this output signal, and then displayed on 
the spectrum analyser, generally in a logarithmic format. 

The output signal is obtained from the following 
equation: 




(8) . 



<p{t) = \27rf(t)dt = —Bl-t 2 



(9) . 



00 



y(t) = x(0 * h used (0 = j h med (r ) • x(t - t) dr 



-00 



Introducing equation (7) leads to the following equation: 
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00 

y(t) = J^(r)V(f-r)«'* ( '- r) </r 

-00 



Introducing equation (9) finally leads to the following 
equation : 

AO = \h med (T)-v\t-T)e * dr 

-co 



The following equation is obtained by multiplying out: 

yit) = ? ,\ h used( T ) eK -v{t-r)e K dr 

v J \ / 

(10) , 

wherein the first term e jlp(t) is not problematic, because 
ultimately |y(t)| will be displayed. The impulse 
response : 



MO = UO-^^' 

(id 



is introduced into the eguation. The index "disp" denotes 
"displayed", because the spectrum of this impulse 
response is displayed, as will be shown below. 
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According to equation (8), the following equation is 
obtained by conversion: 




(12) . 

The following equation is obtained by substitution into 
equation (10) : 

y{t) = J h disp (r) -v\t- t) dr 

-co 

(13) . 

Now, some interesting statements can be made: comparing 
equation (13) with the Fourier analysis in equation (1) 
shows that 

1. The "displayed'"'" resolution fiiLei. udispvw described 
according to equation (11) rather than the used 
resolution filter h used (t) is displayed in the spectrum 
analysis. With a slow sweep for approximately k^2 , h use d(t) 
and h diS p(t) are approximately in agreement. However, with 
a faster sweep, significant differences occur. In this 
case, the level collapses, and the displayed resolution 
filter becomes broader (the filter can no longer build 
up) . 

2. In equation (13), by contrast with the Fourier 
analysis, the time signal displaced by t is used rather 
than v(x). As a result, the spectrum analyser evaluates an 
observation interval sliding over time, which is not in 
itself problematic. However, the question regarding the 
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influence of the rate of the sliding observation window 
on the output spectrum is particularly relevant. 

In order to evaluate the question of the sliding 
observation window in point 2 in an improved manner, it 
is advisable to apply the Parseval's theorem, as set out 
below, to equation (13) : 



oo 



oo 



J x x (t) • x* 2 (r) dr 



J X x {F)>X* 2 (F)dF 



-00 



-00 



By substitution of: 



r jO)(t)T 



4 X X (F) = H di (F + f(t)) 



T 



JIk F t 



x 2 (t) = v(t-r) 



-+ X 2 (F) = V(-F)-e 



equation (13) can be described by the following 
equations : 



00 




-00 



00 




\H disp (F-f(t)).V\F) e^ Ft dF 



-00 



(14) . 



Accordingly, as anticipated, a convolution of the input 
spectrum with the resolution filter is obtained according 
to the following equation: 
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y(t)= e™ H disp {f(t))*[V{f(tj)e-*"^\ 

By inserting equation (12) into equation (14), the 
following equation is finally obtained: 

y(t) = e Mt) \H disp {F-f{t))-V\F)e *- dF 

-oo 

(15) . 

DE 101 05 258 Al derives only a sweep-optimised Gaussian 
filter. This sweep-optimised Gaussian filter must be a 
linear-phase filter . 

New insights have shown that the sweep-optimised filter 
can be derived for a random filter. The given filter can 
be random both with regard to its magnitude and also with 
regard to phase. It is particularly relevant that no 
restrictions are placed on phase. In the case of the 
sweep filter known from DE 101 05 258 Al, this degree of 
freedom could not be exploited, because the linear-phase 
properties were reguired in that context. With a 
randomly-specifiable phase, minimal-phase filters can now 
be realised according to the invention, which are 
optimised with regard to the necessary transient time. 

Eguation (10) defines the "displayed" impulse response 
hdisp(t) and the "used" impulse response h use d(t). In this 
context, the impulse response h d i Sp (t) describes the 
transform of the frequency response H d i Sp (f) of the 
resolution filter displayed in the sweep, while the 
impulse response h use d(t) is the transform of the used 
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filter with the frequency response H use d(f). The 
correlation between the two impulse responses: 

is already known from equation (11) . Now, the procedure 
according to the invention consists in developing the 
frequency response of the displayed filter Hdi S p(f) rather 
than the frequency response H useC i(f) of the used filter. 
Since only the modulus frequency response is displayed in 
spectrum analysis, the phase can be selected at random. 
According to this design, the impulse response h disp (t) is 
calculated by re-transformation. In the next stage, the 
sought impulse, response h use d(t) of the sweep-optimised 
filter is calculated according to the above formula with 
the sweep rate k. 

The individual design stages are described in greater 
detail below: 

1. Specification of the desired displayed modulus 
frequency response |Hdi sp (f) I : 

A Gaussian filter is often used. However, filters with a 
less steeply-declining modulus frequency response can 
also be of interest, because this can reduce the group 
delay response. This is associated with a shorter 
transient time, which is particularly desirable in the 
case of applications with frequent transient processes of 
the filter. 

2. Specification of the phase of H disp (f): 
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In principle, the phase can be specified randomly. In 
order to achieve a minimal group delay response, it is 
advisable to use a minimal-phase filter. Details of the 
design are given below. In this context, it is assumed 
that the filter is realised as a digital filter, that is 
to say, the discrete impulse response h d i S p(k) with 
k= [0, nof Tap g-l] is calculated, wherein nof Ta ps describes the 
number of taps. 

The simplest case is a design of a linear-phase filter 
with a length nof Ta ps with a specified modulus frequency 
response |Hdi sp (f)l- Currently available methods include 
the Remez algorithm and the MMS algorithm. Following 
this, the zero positions of the transmission function in 
the complex z-plane are determined. The zero positions 
outside the circle of radius 1 are then reflected into 
the circle of radius 1, so that the modulus frequency 
response remains unchanged. This method is sub-optimum, 
because it always leads to minimal-phase filters with 
doubled ^ciu positions, which limits the degree of 
freedom. 

Figure 3 shows an example of the discrete impulse 
response h(k) as a function of the sampling index of a 
minimal-phase Gaussian filter by comparison with a 
linear-phase Gaussian filter. It is evident that the 
minimal-phase filter has a substantially shorter group 
delay response than the linear-phase filter. Of course, 
the minimal-phase filter has the same number of taps 
nof Taps =161 as the linear-phase filter, that is to say, 
the transient time is of the same length. However, during 
the sweep, the transient error of the minimal-phase 
filter decays substantially more rapidly than in the case 
of the linear-phase filter, as can be seen from a 
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comparison of Figure 4 with Figure 5. Figure 4 shows the 
sweep with a minimal-phase filter, while Figure 5 shows 
the same sweep with a linear-phase filter. In each case, 
the input signal is a discrete spectral line. It is 
evident that, in the case of the minimal-phase filter, 
the error has already decayed so severely during the 
transient phase, i.e. immediately after the "first lobe", 
that approximately the second half of the transient phase 
can be used for the analysis. By contrast, with the 
linear-phase filter, the analysis cannot be started until 
later . 

An improved filter design is achieved with the following 
method. The linear-phase filter is not designed with 
nof Ta ps taps, but rather with double the length 2*nof Ta ps- 
As a target function, |H disp (f) 1 2 is also specified rather 
than |H disp (f) I . The same applies for any cost functions 
used. Following this, the zero positions in the circle of 
unit radius are calculated from the determined digital 
filter h^\k) . Tne mirror-image symmetrical zero 
positions outside the circle of unit radius are rejected. 
The impulse response h disp (k) generated in this manner has 
the desired number of taps nof Taps . Moreover, it also has 
the desired target modulus frequency response |H disp (f)|. 
The filter calculated in this manner has no doubled zero 
positions, so that this design exploits the full degree 
of freedom. 

3. Calculation of h used (k): 

The impulse response h disp (k) is now available. According 
to equation (11), the impulse response of the filter to 
be used is calculated using the following equation: 
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-j—B 2 (k-k ) 2 -T 2 

KedW) = h di Sp ( k )- e K ^ ° with k = [0,nof Taps - 



(16), 



wherein T a is the sampling period of the digital filter. 
By contrast with equation (11), the parameter k 0 , which 
brings about a displacement of the spectrum within the 
spectral range, has been introduced according to the 
invention. Various desirable effects can be achieved with 
the parameter k 0 : 



1. Compensation of the frequency overshoot: as a result 
of the group delay of the resolution filter, the 
corresponding frequency of the output signal has a 
similar overshoot. This effect can be compensated by 
setting an appropriate k 0 . 



2. Minimisation of the required bandwidth of H usec i(f)- 
with an increasing sweep rate (smaller k) , the middle of 
the spectrum H useC i(f) is no longer at the frequency origin 
f=0, but is displaced towards the higher frequencies. 
With the present digital system, the required sampling 
rate f a would therefore have to be unnecessarily 
increased, in order to continue to fulfil the sampling 
theorem. This effect can be avoided by an appropriate 
choice of k 0 . In a good approximation, k 0 =k max should be 
selected, k max being the time of the maximum of |h disp (k) | . 

Reference must also be made to a procedure, which is 
meaningful for the implementation: h use d(k) should never be 
pre-calculated for different sweep rates k and stored in 
the device; only h disp (k) should be pre-calculated. The 
relevant impulse response h use d(k) will be calculated in 
the device according to the specifications in equation 
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(16) after the desired sweep rate is known. This method 
has the advantage that the housekeeping and memory 
requirement for the individual impulse responses does not 
occur. Furthermore, the impulse response for the present 
k is always used. As a result, no quantisation errors 
occur, because the impulse response h use d(k) for the exact 
k-f actor set is not available. 

The closed display of the impulse response of the complex 
resolution filter is determined as follows. 

For the derivation of the closed display of the impulse 
response h usec i of the filter, the time-continuous case will 
first be considered. In this case, the free variation 
parameter is to=k 0 -T a . 

A Gaussian resolution filter with the bandwidth B res is 
used for the spectrum analysis. The ^displayed" 
resolution filter should have the following impulse 
response ana transmission function: 



The amplitude-corrected display of the spectral lines is 
secured with | H disp ( f =0 ) | =1 . The phase characteristic cp(f) 
of the transmission function can be selected randomly, 
because |Hdi sp (f) I is displayed in the spectrum analysis. 
The phase response <p(f) can be determined, for example, 
according to the method described above, so that the 







(17) 
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transmission function is a minimal-phase function, which 
leads to a rapid transient of h d i SP (t). The Fourier re- 
transform of e j(p(f) is haiip(t) . Since a multiplication with 
e j(p(f> is implemented in the frequency domain, a 
corresponding convolution with the impulse response 
haiip(t) of an all-pass filter must be implemented in the 
time domain, which is symbolised by the character *. 
Because of the correlation indicated in equation (11) 
between h disp (t) and h US ed(t), q>(t) describes not only the 
phase response of the transmission function H disp (f), but 
also the phase response of the transmission function 
Hused(f) of the resolution filter. 



The impulse response of the resolution filter to be used 
is obtained according to equation (16) in a time- 
continuous display as: 

-j—B 2 res {t-t Q f 

MO = h disp (t)-e K 

(18) . 

The following equation is obtained by substitution: 



Ksed (0 -\/ 2In(2) ^ res 



21n(2) 



T 

V J res j 



MO 



(19) 



The transition to the discrete impulse response follows 
from: 



KseAk) = T a h M (t=kr a ). 
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By substitution of equation (19) , the following equation 
is obtained: 



used 



(k) = EEL 

v ' V 2In (2) 



B T 

res a 



— (-1 



-y^B«c-'o) 2 

9 A 



t=kT n 



(20) 



With T res =l/B re s, B res =resolution bandwidth with a 3dB 
signal decay relative to the maximum and T a = sampling 
period in the baseband. 

The following general presentation for the impulse 
response is obtained : 



*u*(*) = c, 



-jC^k-kvf-T? 



In Lais eqaation, k denotes the campling index, and T- 
denotes the sampling period. Ci, C 2 and C 3 are constants, 
wherein the value of the constant Ci is preferably: 



n 



- 1 = \|21n(2) 



and wherein B res is the bandwidth of the resolution 
filter. 

The value of the constant C 2 is preferably: 



C 2 = 



1 

21n(2)V 



2 , 
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wherein T res =l/B res is the reciprocal bandwidth B res of the 
resolution filter. 



The value of the constant C3 is preferably 



£. " 2 
K 



wherein B res is the bandwidth of the resolution filter, 
and K is the k-factor of the resolution filter, which is 
defined by the equation: 



and f (t) is a linear frequency variable with time t, 
which is supplied to the mixer 3 of the spectrum analyser 
connected upstream of the resolution filter 4. 



However, the constants Ci, C 2 and C 3 may also conceivably 
be specified in a different manner within the framework 
of the present invention. 

While the present invention has been described in 
connection with a number of embodiments and 
implementations, the present invention is not so limited 
but covers various obvious modifications and equivalent 
arrangements, which fall within the purview of the 
appended claims. 



